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Abstract 



We report on global fits of optical-model parameters to 90 data points for 
p X-rays and 17 data points of radiochemical data put together. By doing 
separate fits to the two kinds of data it is possible to determine phenomeno- 
logically the radial region where the absorption of antiprotons takes place and 
to obtain neutron densities which represent the average behaviour over the 
periodic table. A finite-range attractive and absorptive p-nuclear isoscalar 
potential fits the data well. Self-consistent dynamical calculations within the 
RMF model demonstrate that the polarization of the nucleus by the atomic 
antiproton is negligible. 
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I. INTRODUCTION 



Experiments on strong-interaction effects in antiprotonic atoms usually involve measure- 
ments of X-rays from p-atomic transitions, providing values of level shifts and widths, the 
latter either measured directly or indirectly via the observed yields. Global analyses of large 
data sets, in terms of an optical model approach, have provided parameters for p-nucleus op- 
tical potentials [1,2] which fit well all the data then available, close to 50 data points, across 
the periodic table. It was found that the radial region which is 'sampled' by the p-atom data 
is approximately between the half-density nuclear radius and some 4 fm outside of it. The 
data could be well fit by a strongly attractive and strongly absorptive isoscalar p-nuclear 
potential V opt , underlined by a NN s-wave interaction, in which the imaginary (absorptive) 
part outweighs the real (attractive) part. Gradient (p-wave) terms in V op t were not required 
by the p-atomic data. Furthermore it was concluded, for any global fit of p-atom data, that 
the isovector components of V opt were relatively insignificant, provided that single-particle 
densities which describe correctly the nuclear surface region and beyond were used. We 
note that since V opt is strongly absorptive, its extrapolation into nuclear-matter densities 
is highly model dependent. Even if a unique extrapolation were possible, the existence of 
well-defined nuclear bound states of antiprotons would have been unlikely due to the large 
widths expected. 

The experimental situation has changed significantly with the availability of high-quality 
data for several sequences of isotopes along the periodic table due to the PS209 collaboration 
[3] which made it possible to perform global fits to larger and more accurate data bases 
consisting of close to 100 data points. A good description of these new data was achieved 
in global analyses by using a finite-range V opt in which a pN s-wave interaction is folded 
with the nuclear density p = p n + p p , as reported preliminarily in Ref. [4]. Another type of 
experimental information on the p-nucleus interaction in p atoms is obtained using the so- 
called radiochemical method. By studying the production of nuclei differing from the target 
nucleus by the removal of one neutron or one proton it is possible to obtain information on 
the ratios of p absorption on neutrons to p absorption on protons in the periphery of the 
target nucleus, at a distance about 2.5 fm beyond the half-density charge radius according 
to Ref. [5]. At such distances the nuclear density is dominated by the neutron density 
which is not directly measured or model-independently determined, and it is vital to use 
reliable nuclear models rooted in a sound nuclear phenomenology to extract these densities 
for use in optical-model calculations. Close to 20 target nuclei have been studied by the 
radiochemical method [6,7] and in several cases [8-10] good consistency was found between 
the radiochemical data and the X-ray data, within an optical-model approach. 

The present paper reports on global fits of optical potentials to all those experimental re- 
sults, consisting of 90 points of X-ray data and 17 points of radiochemical data put together. 
Our aim in this paper is to determine the p-nuclear interaction potential at the nuclear sur- 
face and beyond, while identifying simultaneously the largely unknown neutron densities 
which provide a good fit to the combined data in terms of V opt that depends sensitively on 
these densities. The organization of the paper is as follows. In section II we verify that 
the conventional approach of static atomic calculations, i.e. using static nuclear densities, 
is valid. This is achieved by performing dynamical calculations within a relativistic mean 
field (RMF) approach where the nucleus is allowed to be polarized by the atomic p. Note 
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that for p-nuclear configurations, in contrast, dynamical RMF calculations for antiprotons 
[11] have yielded considerable nuclear polarization effects. Section III presents the p-nucleus 
potential and the methodology of the present work, in particular as regarding the use of 
neutron densities. In section IV we present results, first only for the X-ray data, followed 
by results for the radiochemical data where we examine under what conditions the two sets 
of results are consistent with each other. Global fits to the combined data base of X-ray 
and radiochemical data are then presented and conclusions are drawn on the admissible 
neutron-density shapes. Section V summarizes the findings of this work. 



II. DYNAMICAL CALCULATIONS OF ANTIPROTONIC ATOMS 

Self-consistent RMF calculations of deeply bound p - nuclear states [11], as well as K- 
nuclear states [12], revealed extremely large dynamical polarization effects. The core nucleus 
is highly compressed and its binding energy significantly increases due to the strongly attrac- 
tive p - nucleus interaction. It is therefore mandatory to check to what extent polarization 
effects are important in the calculations of p- atomic states and whether the static approach 
hitherto used in atomic calculations is justified. In order to address this problem, we per- 
formed self-consistent dynamical calculations of p - atomic states within the RMF model. 

In the RMF framework the core nucleons and the p interact through the exchange of 
isoscalar-scalar (<r) and -vector (u), isovector- vector (p), and Coulomb (A) fields, which 
are treated in the mean-field approximation. The advantage of the RMF approach is that 
antinucleons (N) are naturally included in the Lagrangian density and, consequently, the 
derivation of the relevant equations of motion is straightforward. The Dirac equation for 
nucleons and antinucleons (i = N, N) reads: 

-id- V + P(Mi + Si) + K] ^ a = e-^ a , (1) 

where Si = g a iO~, Vi = + g P i T 3P\ + Ci^-^-A , and a denotes quantum numbers of 

single-particle states. 

The presence of the p in the nuclear system under consideration modifies the source 
terms in the equations of motion for the meson fields: 

(-A + ml + da + c 2 o- 2 ) a = -g cN psN - QgnPsn , 

(-A + ml + duj 2 ) LO = QuNpBN + QuiNpBN , (2) 

(-A + mfj pi = g pN p 3N + g P Np 3 N , 
(-A) A = e p p p + e p p p , 

where psi, PBi, P3i, Pj (j = P, P) are the scalar, vector (baryon), isovector, and proton and 
antiproton densities, respectively. 

While for nucleons the Dirac equation (1) was solved, for the antiproton we actually 
solved a Schrodinger equation with a complex p optical potential. The real part of the p 
potential, ReV opt , was constructed as the Schrodinger equivalent potential from the scalar 
and vector meson mean fields: 

ReK P t = S R + V R + (S% - V£)/2M . (3) 
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Note that, for antiprotons, both Sfj and are attractive owing to G-parity. When solving 
Eqs.(2) for the meson mean fields, we make the approximation p s ^- = p B fj = p^f, where p^ is 
the p density formed by solving the Schrodinger equation (3). This approximation amounts 
to neglecting a few percent difference at most between the scalar and vector densities, a 
neglect that turns out to be completely insignificant for the dynamical effects calculated in 
p atoms. The p optical potential is dominated by its imaginary term [2] which is due to 
the annihilation of antiprotons on nucleons. Im\4 pt was taken proportional to the nuclear 
density p(r) and its depth was adjusted to fit the atomic data. It is to be noted that in 
this dynamical approach the nuclear density p is not a static quantity but is affected by the 
antiproton interacting strongly with the core nucleons via meson mean fields. 

In our dynamical calculations of the p-nucleus systems, the standard RMF Lagrangian 
with the nonlinear parameterizations TM1 [13] and due to Sharma et al. [14] were used. 
Applying the G-parity transformation as a starting point to determine the antinucleon (p) 
couplings gjff (j = a, u, p) but realizing that the presence of strong annihilation channels 
and various many-body effects could cause significant deviations from the G-parity values, 
we introduced following Ref. [11] a reduction parameter e (e < 1) to allow for departures 
from the ideal G-parity limit: 

9aN — £9aN, 9ujN = —^9ujN , 9pN = ^9pN ■ (4) 

The reduction parameter e was determined by fitting the calculated energy shift and width 
of a particular p atomic state to the corresponding experimental result. Static calculations 
showed that fits to the atomic p data for 16 0, 40 Ca and 208 Pb produced very nearly the same 
potential parameters as obtained from fits to the full data base (but with larger errors with 
which the optical-potential parameters are determined). Consequently we chose these three 
nuclei for checking polarization effects by studying the 3d atomic state of 16 0, the 4/ state 
of 40 Ca and the 9k state of 208 Pb. The value of the reduction parameter e needed to fit the 
atomic data was between 0.15 ( 40 Ca) and 0.35 ( 208 Pb). This indicates a large deviation of 
the p couplings from the G-parity values. 

The coupled system of equations (1) for nucleons, (2) for the meson- and electromagnetic 
fields, and the Schrodinger equation for p was solved self-consistently by iterations. In 
order to check the numerical procedure, p-nuclear states were calculated using the TM1 
parameterization and the results of Ref. [11] were reproduced. Large polarization effects 
were found for p-nuclear states, as expected. Subsequently, dynamical effects for p- atomic 
states were studied. Comparing the energy shifts and widths calculated self-consistently 
with the predictions of the corresponding static calculations, we found totally negligible 
effects of less than 1 eV for the above three p-atomic states. 

Finally, deeply bound p-atomic states [15] were also calculated in 16 0, 40 Ca, and 208 Pb. 
For illustration, Table I presents calculated binding energies (B p ) and widths (T p ) for the 
Is and 2p p-atomic states in Ca, comparing between the results of dynamical and of static 
calculations. The differences of few eV, caused by the polarization effects, are exceedingly 
small and experimentally unobservable. We conclude that static analyses of p-atomic data 
which neglect the rearrangement of the core nucleons are adequate. 
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III. METHODOLOGY 



In order to preserve the connection to previous studies of hadronic atoms [2] , the Klein- 
Gordon (KG) equation is used in the form 



Here, V c denotes the static Coulomb potential for the p due to the finite charge distribution 
of the nucleus, including the first-order a (Zen) vacuum-polarization potential, ji is the p- 
nucleus reduced mass and B = Bp + iTp/2 is the complex binding energy. The interaction 
of antiprotons with the nucleus is described here in terms of an optical potential V opt which 
in the simplest 'tp' form is given by 



where p n and p p are the neutron and proton density distributions normalized to the number 
of neutrons N and number of protons Z, respectively, A = N + Z, and M is the mass 
of the nucleon. In the zero-range 'tp' approach the parameters &o an d &i are minus the p- 
nucleon isoscalar and isovector scattering lengths, respectively, otherwise these parameters 
may be regarded as 'effective' and are obtained from fits to the data, either as given above 
or in a finite-range folding model [4]. The explicit inclusion of an isovector term in V opt 
is motivated by the fact that much data now exist [3] for chains of isotopes and because 
the radiochemical method provides information on ratios of neutron to proton densities. A 
least-squares search procedure was used to adjust the potential parameters so as to obtain 
a best fit to the experimental shift and width measurements. 

A note on wave equations is in order here. The quadratic term (V c + B) 2 in the KG 
equation (5) has little effect on the calculated strong-interaction shift and width, and omit- 
ting it gives rise to the nonrelativistic Schrodinger equation. One might ask why use the 
KG equation as extension into the relativistic domain instead of the Dirac equation for the 
p fermion. Indeed when interpreting experimental transition energies in order to extract 
the strong interaction effects it is essential to use the Dirac equation with finite size nuclear 
charge distribution and vacuum polarization terms [16]. However, strong interaction effects 
in p atoms are normally given as the proper averages over the fine structure components. 
The use of the KG equation rather than the Dirac equation is numerically justified when 
fine-structure effects are negligible or are treated in an average way, as for the X-ray transi- 
tions considered here. The leading j dependence (j = /± |) of the energy for solutions of the 
Dirac equation for a point-charge 1/r potential goes as (j + and on averaging it over 
the projections of j gives rise to (/ + which is precisely the leading / dependence of the 
energy for solutions of the KG equation. The higher-order contributions to the spin-orbit 
splitting are suppressed by 0(Za/n) 2 which is of order 1% for the high-n X-ray transitions 
encountered for antiprotons. We checked numerically for few typical cases that the spin-orbit 
averaged shifts and widths thus obtained differ by less than 1% from the (2j + l)-average 
of the corresponding quantities obtained by solving the Dirac equation. This difference is 
considerably smaller than the experimental errors placed on the measured X-ray transition 
energies and widths. 



A - 2fi{B + K P t + V c ) + (V c + Bf i> = (h = c = 1). 



(5) 




(j, A - 1 
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The nuclear densities are an essential ingredient of the optical potential. The density 
distribution of the protons is usually considered known as it is obtained from the nuclear 
charge distribution [17] by unfolding the finite size of the charge of the proton. The neutron 
distributions are, however, generally not known to sufficient accuracy. A host of different 
methods have been applied to the extraction of rms radii of neutron distributions in nuclei 
but the results are sometimes conflicting, see Refs. [18-23]. For many nuclei there is no 
direct experimental information whatsoever on neutron densities and one must then rely 
on models. To complicate things further we note that there is a long history of conflict 
between values of neutron rms radii derived from experiments using hadronic projectiles 
and neutron rms radii obtained from theoretical calculations. For that reason we have 
adopted a semi-phenomenological approach that covers a broad range of possible neutron 
density distributions. 

Experience with pionic atoms showed [24] that the feature of neutron density distribu- 
tions which is most relevant in determining strong interaction effects in pionic atoms is the 
radial extent, as represented for example by r n , the neutron density rms radius. Other fea- 
tures such as the detailed shape of the distribution have only minor effect. For that reason 
we chose the rms radius as the prime parameter in the present study. Since r p , the rms 
radius for the proton density distribution, is considered to be known, we focus attention on 
values of the difference r n — r p . In previous analyses of p radiochemical data [8,23] a linear 
dependence of r n — r p on (N — Z) / A was employed, namely, 

N ~ z 

r n -r p = a— - — + 7 , (7) 

with a close to 1.0 fm and 7 close to zero. The same expression with a close to 1.5 fm 
was found [24] to represent well results of RMF calculations [25] for stable nuclei, but these 
values of r n — r p are larger by about 0.05 — 0.10 fm than the 'experimental' values in medium- 
weight and heavy nuclei used in recent relativistic Hartree-Bogoliubov (RHB) versions of 
mean-field calculations [26,27]. Expression (7) has been adopted in the present work and, 
for lack of better global information about neutron densities, the value of a was varied over 
a reasonable range in fitting to the data. This procedure is based on the expectation that 
for a large data set over the whole of the periodic table some local variations will cancel out 
and that an average behaviour may be established. 

In order to allow for possible differences in the shape of the neutron distribution, the 
'skin' and the 'halo' forms of Ref. [8] were used, as well as an average between the two. 
Adopting a two-parameter Fermi distribution both for the proton (unfolded from the charge 
distribution) and for the neutron density distributions 

PnAr) ~ l + eM(r-R n , P )/a n , p ) ' W 

then for each value of r n — r p in the 'skin' form the same diffuseness parameter for the 
protons and the neutrons a n = a p was used and the R n parameter was determined from the 
rms radius r n . In the 'halo' form the same radius parameter R n = R p was assumed and 
was determined from r n . In the 'average' option the diffuseness parameter was set to be the 
average of the above two diffuseness parameters a^ vc = (a p + a^)/2 and the radius parameter 
R n was then determined from the rms radius r n . In this way we have used three shapes of 
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the neutron distribution for each value of its rms radius all along the periodic table. Figure 
1 shows as an example the densities for 120 Sn after a finite range folding (see below). The 
relevant radial region for the p-nucleus interaction is 8-9 fm in this case, and it is seen that 
the neutron densities are quite different for the three models. 

Returning to the optical potential Eq.(6), we note that in Ref. [4] it was found that the 
X-ray data could be described better with a finite-range interaction, where the densities of 
Eq.(6) were replaced by 'folded' densities 

/(r) = / dr'p^J^e-^lf , (9) 

with the value of (5 — 0.85 fm used for both the real and imaginary parts of the potential. 
For the analysis of the radiochemical data we adopt the approach of Refs. [6-8] that the 
method is sensitive to the neutron to proton density ratio close to 2.5 fm outside of the half- 
density radius of the charge density [5]. The experimental ratios of absorption on neutrons 
to absorption on protons were therefore compared to 

Im(6 + 6i)/ w 
Im(6 - bi)I p 

where I niP are the volume integrals of the neutron and proton densities, respectively, either 
between 2.0 and 3.0 fm or between 2.5 and 3.5 fm outside of the half-density radius of the 
charge density. For a finite-range potential the folded densities were used. No use was 
made of atomic wavefunctions in calculating the ratios because their effect largely cancel 
out in the ratios. Choosing the range of integration was guided by the conclusions of Ref. 
[5] which obviously were based on atomic wavefunctions. By requiring consistency between 
radiochemical data and X-ray data we could independently check the conclusions of Ref. [5] 
as explained below. 
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IV. RESULTS AND DISCUSSION 
A. X-ray data 

A brief report on global fits to the present extended X-ray data has already been pub- 
lished [4]. The data base included strong interaction level shifts and widths for antiprotonic 

atom levels in 16 < 18 0, 40,42,43,44,48 Ca> 54,56,57,58^ 58,60,62,64^ 90,96^ 106,116^ 112,116,120,124^ 

i22,i24,i26,i28,i30 Te and 2os pb; a total of gg points [3]. i t was f omK } t h at by introducing a finite 
range into the interaction within the folding-model approach of Eq.(9) consistently better 
fits to the data were obtained compared to the corresponding zero-range approach. 

Figure 2 shows values of x 2 obtained for the three shapes of neutron distributions dis- 
cussed above as function of the slope parameter a in the expression for r n — r p , Eq.(7), using 
a finite-range parameter of (3=0.85 fm. It is seen that both the 'halo' model and the 'aver- 
age' model for p n produce very acceptable fits with x 2 per point of about 2.2. Recall that 
values of a are likely to be between 0.9 and 1.3 fm, as the value of 1.5 fm found from fits to 
results of RMF calculations is now regarded as an over-estimate. The results for the 'halo' 
version confirm earlier findings [8] that for this shape of neutron density agreement with 
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experiment is obtained for a very close to 1 fm. For the 'skin' variety of p n the minimum 
of x 2 is obtained for a value of a close to 2 fm which is unreasonably large [23,26,27]. The 
'average' type of p n is quite acceptable both in the quality of fit and in the value of a at the 
minimum of x 2 ■ For this model the slope of p n is smaller at large radii than the slope of p p , as 
expected on grounds of binding energies, yet the model is not restricted in having the same 
half-density radii for the proton and for the neutron distributions. Note that from global fits 
to strong interaction shifts and widths in pionic atoms it is found [24] that the 'halo' shape 
is rejected, whereas the other two shapes are quite acceptable. However, the pionic atom 
data are sensitive to considerably smaller radii (and consequently larger densities) than is 
the case for p atoms. Figure 3 shows values of the isoscalar amplitude b from fits to the 
data for the 'average' type of neutron density using a finite-range interaction. The isovector 
amplitude was found to be consistent with zero [4] and we return to this point below. 



B. Radiochemical data 

Experimental ratios of p absorption on neutrons to absorption on protons were taken 
from Refs. [6,7]. Initial calculations showed that very large contributions to the resulting 
X 2 came from 106 Cd and 112 Sn and subsequently these two nuclei were excluded from the 
data set. Possible explanations for the problem with these two nuclei in terms of a pp quasi 
bound state are given in Ref. [28]. We have therefore used 17 values of absorption ratios 
for the following nuclei: 48 Ca, 58 Ni, 96 Zr, 100 Mo, 96 < 104 Ru, 116 Cd, 124 Sn, 12 8>i30 Te) i44,i54 Sm> 
148 Nd, 160 Gd, 176 Yb, 232 Th and 238 U. With the potential parameter &i consistent with zero 
the ratios Eq.(lO) become independent of the parameters of the potential, but they are 
quite sensitive to values of r n — r p . This is seen clearly in Fig. 4 where x 2 values for the 
radiochemical data are shown as function of a for the 'halo' and for the 'average' types of 
neutron densities. With the proton density held fixed the calculated ratios obviously depend 
on the neutron density and on the radius where the ratios of densities are calculated. To 
avoid the use of densities at a point we have integrated the densities over a distance of 1 fm, 
and Fig. 4 shows results for integrations between 2.5 and 3.5 fm and between 2 and 3 fm 
beyond the half-density radius of the charge distribution. However, due to the exponential 
decrease of the densities at such large radii the integrals are dominated by the densities close 
to the lower limit of the range of integration. 

Comparing positions of the minima in Fig. 4 for the radiochemical data with Fig. 2 
for the X-ray data, we note a consistency between the two regarding the parameter a when 
the integration range for the radiochemical data is 2.5 to 3.5 fm beyond the charge radius, 
whereas moving that range down by 0.5 fm leads to inconsistency between the two types 
of data. This is a significant test because both types of data depend differently on the 
nuclear densities. This result confirms in a phenomenological way the theoretical conclusion 
of Wycech et al. [5] that most of the absorption takes place close to 2.5 fm beyond the charge 
radius. We therefore adopt the 2.5 to 3.5 fm segment as the integration range for the global 
analysis of the combined X-ray and the radiochemical data. 
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C. Global fits to X-ray and radiochemical data 

Fits to all the 107 data points due to the X-ray and radiochemical data were made, using 
the various shapes for the neutron densities as described above and varying the isoscalar 
potential parameter b and the isovector parameter b±. Finite-range folding was assumed 
throughout with a range parameter (3=0.85 fm. For the radiochemical data the integration 
range was chosen between 2.5 and 3.5 fm beyond the charge half-density radius, as described 
above. The real part of b\ was always found to be consistent with zero and subsequently 
it was excluded from the fits. Figure 5 shows results for the 'average' shape of the neutron 
densities and Fig. 6 shows similar results for the 'halo' shape. In the lower part of the 
figures we show values of y 2 for the X-ray data separately and for the combined X-ray and 
radiochemical data. The upper part shows the potential parameters. 

It is seen that although Im&i assumes non-vanishing values away from the minima of 
X 2 , at the minima this parameter is essentially zero. Accepting those minima for the 'halo' 
and for the 'average' shapes with a between 0.9 and 1.3 fm as a fair representation, on the 
average, of neutron densities, we end up with the values of b as found in Ref. [4], namely 
Re&o=l-3±0.1 fm, Im&o=l-9±0.1 fm, with a folding range of (3=0.85 fm. 

The vanishing of Im&i at the best-fit points for the two acceptable models for neutron 
densities deserves a comment. We note that the four NN potentials considered in Table 3 
of Ref. [1] yield for Im&x values which are a factor of 10-20 smaller than the corresponding 
values of Im6 , so it is conceivable that also for the effective phenomenological parameters 
a similar situation will hold. With the many groups of isotopes in the present data base 
the dependence on iV — Z must play a role. The neutron densities are obviously properly 
normalized, but the rms radii are not set a priori and they are found, on the average, in the 
process of \ 2 fits. The vanishing of \mbi is presumably an indication of the general validity 
of the neutron densities used here. As a test we note that for the unacceptable 'skin' model, 
where a is close to 2 fm at the minimum, a non- vanishing value is found for Im&i, which 
presumably compensates for the unphysical neutron densities. 

The fits shown in Figs. 2-6 can be marginally improved, without changing the conclu- 
sions of the present work, if the parameterization of r n — r p in Eq. (7) is slightly modified to 
incorporate the expected A dependence of the parameter 7 due to the increased Coulomb 
repulsion within the N = Z core. To this end we have used 7 = —0.0162 A 1 ^ 3 fm to inter- 
polate between r n — r p values in 16 0, 40 Ca and the Sn and Pb isotopes within recent RHB 
calculations [26] which yield values for the (N — Z)/A slope parameter a ~ 1.5 — 1.6 fm. 
This analysis, again, strongly suggests that the prefered neutron densities are best described 
by the 'average' type shape. A more expanded discussion of these results, plus related ones 
for pionic atoms, are relegated to a separate publication. 

V. CONCLUSIONS 

In conclusion, self-consistent dynamical calculations within the RMF model demonstrate 
that the polarization of the nucleus by the atomic antiproton is negligible, thus confirming 
the validity of the common static approach in which the nuclear polarization is disregarded. 
We have performed, for the first time, global fits of optical model parameters to 90 data 
points for antiprotonic X-rays and 17 data points of radiochemical data put together. With 
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the help of separate fits to the two kinds of data we could verify in a phenomenological 
way the theoretical prediction [5] for the radial region where most of the absorption of 
atomic antiprotons takes place. Despite the use of many groups of isotopes all along the 
periodic table, the parameters of the isovector part of the potential could not be determined 
in global fits, and are therefore assumed to be consistent with zero. The optical potentials 
are determined at radii where the density is 2-3% of the central nuclear density. It is 
therefore rather meaningless to extrapolate these potentials to the nuclear interior. However, 
if extrapolated with the nuclear densities used in the present work, then at the center of 
the nucleus the potential is attractive with about 110 MeV depth, and its absorptive part is 
close to 160 MeV deep. While this attraction is fairly strong, about twice that for nucleons 
in nuclear- matter densities, it is much weaker than naive G-parity arguments suggest. The 
world's data on p-nucleus interaction potential are well accounted for with the effective 
parameters Re6o=1.3±0.1 fm, Im6 =1.9±0.1 fm, with a folding range of /?=0.85 fm. 

This work was supported in part by the Israel Science Foundation grant 131/01 and by 
the GA AVCR grant IAA1048305. 



9 



REFERENCES 



[1] C.J. Batty, E. Friedman, A. Gal, Nucl. Phys. A 592 (1995) 487. 
[2] C.J. Batty, E. Friedman, A. Gal, Phys. Rep. 287 (1997) 385. 

[3] A. Trzcihska, J. Jastrz§bski, T. Czosnyka, T. von Egidy, K. Gulda, F.J. Hartmann, 
J. Iwanicki, B. Ketzer, M. Kisielinski, B. Klos, W. Kurcewicz, P. Lubihski, P.J. Napi- 
orkowski, L. Piehkowski, R. Schmidt, E. Widmann, Nucl. Phys. A 692 (2001) 176c. 

[4] E. Friedman, A. Gal, Nucl. Inst. Methods in Phys. Research B 214 (2004) 160. 

[5] S. Wycech, J. Skalski, R. Smolanczuk, J. Dobaczewski, J. Rook, Phys. Rev. C 54 (1996) 
1832. 

[6] P. Lubihski, J. Jastrzebski, A. Trzcihska, W. Kurcewicz, F.J. Hartmann, R. Schmidt, 

T. von Egidy, R. Smolanczuk, S. Wycech, Phys. Rev. C 57 (1998) 2962. 
[7] R. Schmidt, F.J. Hartmann, B. Ketzer, T. von Egidy, T. Czosnyka, J. Jastrzebski, 

M. Kisielinski, P. Lubihski, P. Napiorkowski, L. Piehkowski, A. Trzcihska, B. Klos, R. 

Smolanczuk, S. Wycech, W. Poschl, K. Gulda, W. Kurcewicz, E. Widmann, Phys. Rev. 

C 60 (1999) 054309. 

[8] A. Trzcihska, J. Jastrzebski, P. Lubihski, F.J. Hartmann, R. Schmidt, T. von Egidy, B. 
Klos, Phys. Rev. Lett. 87 (2001) 082501. 

[9] R. Schmidt, A. Trzcihska, T. Czosnyka, T. von Egidy, K. Gulda, F.J. Hartmann, J. 
Jastrzebski, B. Ketzer, M. Kisielinski, B. Klos, W. Kurcewicz, P. Lubihski, P. Napi- 
orkowski, L. Piehkowski, R. Smolanczuk, E. Widmann, S. Wycech, Phys. Rev. C 67 
(2003) 044308. 

[10] B. Klos, S. Wycech, A. Trzcihska, J. Jastrzebski, T. Czosnyka, M. Kisielinski, P. 
Lubihski, P. Napiorkowski, L. Piehkowski, F.J. Hartmann, B. Ketzer, R. Schmidt, T. 
von Egidy, J. Cugnon, K. Gulda, W. Kurcewicz, E. Widmann, Phys. Rev. C 69 (2004) 
044311. 

[11] I.N. Mishustin, L.M. Satarov, T.J. Burvenich, H. Stocker, W. Greiner, Phys. Rev. C 71 
(2005) 035201. 

[12] J. Mares, E. Friedman, A. Gal, Phys. Lett. B 606 (2005) 295. 

[13] Y. Sugahara, H. Toki, Nucl. Phys. A 579 (1994) 557. 

[14] M.M. Sharma, M.A. Nagarajan, P. Ring, Phys. Lett. B 312 (1993) 377. 

[15] E. Friedman, A. Gal, Nucl. Phys. A 658 (1999) 345. 

[16] R. Schmidt, F.J. Hartmann, T. von Egidy, T. Czosnyka, J. Iwanicki, J. Jastrzebski, 
M. Kisielinski, P. Lubihski, P. Napiorkowski, L. Piehkowski, A. Trzcihska, J. Kulpa, R. 
Smolanczuk, S. Wycech, B. Klos, K. Gulda, W. Kurcewicz, E. Widmann, Phys. Rev. C 
58 (1998) 3195. 

[17] G. Fricke, C. Bernhardt, K. Heilig, L.A. Schaller, L. Schellenberg, E.B. Shera, C.W. De 

Jager, At. Data Nucl. Data Tables 60 (1995) 177. 
[18] C.J. Batty, E. Friedman, H.J. Gils, H. Rebel, Advances in Nuclear Physics 19 (1989) I. 
[19] C. Garcia-Recio, J. Nieves, E. Oset, Nucl. Phys. A 547 (1992) 473. 
[20] V.E. Starodubsky, N.M. Hintz, Phys. Rev. C 49 (1994) 2118. 

[21] A. Krasznahorkay, M. Fujiwara, P. van Aarle, H. Akimune, I. Daito, H. Fujimura, Y. 

Fujita, M.N. Harakeh, T. Inomata, J. Janecke, S. Nakayama, A. Tamii, M. Tanaka, H. 

Toyokawa, W. Uijen, M. Yosoi, Phys. Rev. Lett. 82 (1999) 3216. 
[22] B.C. Clark, L.J. Kerr, S. Hamma, Phys. Rev. C 67 (2003) 054605. 



10 



[23] J. Jastrz§bski, A. Trzcihska, P. Lubinski, B. Klos, F.J. Hartmann, T. von Egidy, S. 

Wycech, int. J. Mod. Phys. E 13 (2004) 343. 
[24] E. Friedman, A. Gal, Nucl. Phys. A 724 (2003) 143. 

[25] G.A. Lalazissis, S. Raman, P. Ring, At. Data Nucl. Data Tables 71 (1999) 1. 
[26] T. Niksic, D. Vretenar, P. Finelli, P. Ring, Phys. Rev. C 66 (2002) 024306. 
[27] G.A. Lalazissis, T. Niksic, D. Vretenar, P. Ring, Phys. Rev. C 71 (2005) 024312. 
[28] S. Wycech, Nucl. Phys. A 692 (2001) 29c. 



11 



TABLES 

TABLE I. Binding energies Bp and widths Tp (in keV) of deeply-bound p atomic states in 40 Ca 
calculated within the static and dynamical approaches using the TM1 RMF parameterization. 



nl 


static Bp 


dynamical Bp 


static Tp 


dynamical Tp 


Is 


1835.132 


1835.129 


584.970 


584.967 


2p 


1516.105 


1516.101 


353.299 


353.300 
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FIGURES 




r(fm) 

FIG. 1. Proton and neutron densities for 120 Sn after a finite range folding with j3 =0.85 fm, see 
Eq.(9). Neutron densities are calculated with a = 1.2 fm, see Eq.(7). 
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150 







a (fm) 



FIG. 2. x 2 values from fits to p X-ray data as function of the slope parameter of r n — r p (Eq.(7)) 
for various types of neutron density. 
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12 3 

a (fm) 

FIG. 4. x 2 values from fits to radiochemical data as function of the slope parameter of r n — r p 
(Eq.(7)) for 'halo' and 'average' types of neutron density and for two integration ranges: A for 2.5 
to 3.5 fm beyond the charge radius; B for 2.0 to 3.0 beyond the charge radius, see text. 
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